Study on the operating characteristics of Stirling engine based on an optimized analysis method* 
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The Stirling engine, as a closed-cycle power machine, possesses excellent emission characteristics and broad 
energy adaptability. Second-order analysis methods find extensive use during the foundational design and ther- 
modynamic examination of Stirling engines due to their commendable model precision and remarkable effi- 
ciency. To scrutinize the effect of Stirling engine design parameters on the cyclical work output and efficiency, 
this paper formulates a series of differential equations for the Stirling cycle employing second-order analysis 
methods, subsequently augmenting the predictive accuracy by integrating considerations of loss mechanisms. 
Additionally, an iterative method for the convergence of average pressure is introduced. The predictive capabil- 
ity of the established model is validated via the GPU-3 and RE-1000 experiment data. Per the model, parameters 
such as the operational fluid, the porosity of the regenerator, and the diameter of the wire mesh, and their influ- 
ence on the resulting work output and cyclic efficiency of the Stirling engine are analyzed, thereby facilitating a 
broader understanding of the engine’s functional characteristics. Findings suggest that hydrogen, courtesy of its 
lower dynamic viscosity coefficient, can engender superior output power. The loss due to flow resistance tends 
to increase alongside the rotational speed. Additionally, under conditions of elevated rotational speed, the loss 
from flow resistance declines in cases of increased porosity, and enhancement of the porosity to diminish flow 
resistance losses can boost both the engine’s output work and cyclic efficiency. As porosity rises further, the 
hydraulic diameter and the dead volume in the regenerator continue to expand, causing the pressure drop within 
the engine to become the dominant factor in the output power’s gradual reduction. Furthermore, extending the 
length of the regenerator results in a decrease in output work — though the thermal cycle efficiency initially 
experiences an increase before eventually declining. Based on these insights, the paper pursues optimal designs 


for the Stirling engine. 


I. INTRODUCTION 


As the call for advanced extraterrestrial exploration voy- 
ages perpetually escalates, the confinement posed by solar 
and chemical energy sources becomes glaringly evident in 
cosmic frontier exploration. Nuclear energy possesses the ap- 
titude to furnish sustained power supply, extending from the 
lunar to even the expanse of remote celestial bodies, thereby 
catering to the diverse requisites of space missions, includ- 
ing but not limited to the establishment of enduring outposts 
(Fig.1). Graced with an extended lifespan, exceptional energy 
quantum, and steadfast operational efficacy, the Space Nu- 
clear Reactor Power System (SNRPS)[1] triumphs as the op- 
timal selection for executing interstellar missions. The power 
prerequisites encompassing various tasks[2, 3] of cosmic ex- 
ploration, such as deep space voyages, orbital space stations, 
and planetary surface infrastructures, within the range of 10 
to 1000 kilowatts[4, 5]. For such space-bound expeditions, 
a power conversion system, imbued with a lengthy service- 
able life, high dependability, and substantial power density, 
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is pivotal. The heat pipe-cooled reactor[6—11], an innova- 
tive nuclear apparatus, harnesses heat pipes to transfer the re- 
actor’s core heat directly to its secondary system. The Stir- 
ling engine, crafted by the ingenious British engineer Robert 
Stirling in 1816[12], heralds a plethora of benefits, encom- 
passing adaptability with virtually all energy sources, dimin- 
ished noise levels, and an indifference to pressure fluctua- 
tions. Recent years have bore witness to notable headway 
in the fabrication of kilowatt-scale SNRPS that amalgamates 
Stirling engines with heat pipe-cooled reactors[ 13-16]. The 
Stirling cycle can achieve a relatively high thermodynamic- 
to-electrical conversion efficiency, theoretically ranging from 
20% to 40%[17]. Hence, the integration of heat pipe-cooled 
reactors with Stirling engines has been actively pursued as an 
area of research. The U.S. has harnessed the principles of the 
Stirling cycle in the conceptualization and design of a diverse 
array of kilowatt-scale heat pipe fission reactors, such as the 
Heat Pipe Mars Exploration Reactor (HOMER)[18] and the 
Kilowatt Reactor Utilizing Stirling Technology (KRUSTY, 
(Fig.2))[19-21]. In March 2018, KRUSTY successfully op- 
erated as a fission power system and was the first nuclear- 
powered operation of any truly new reactor concept in the 
United States in over 40 years[22]. The performance of 
the Stirling engine influences the thermal output of the re- 
actor, so it is necessary to analyze its operational character- 
istics. The Stirling cycle thermodynamic model is widely 
used to describe and predict the process and thermal perfor- 
mance of the Stirling cycle. Ideally, the Stirling cycle con- 
sists of two isothermal processes and two constant-capacity 
processes, and the thermal conversion efficiency is equal to 
the Carnot cycle efficiency. Nonetheless, in practical terms, 


Fig. 1. Space nuclear energy 


various energy dissipations demand consideration, an aspect 
of paramount importance in the construction of an accurate 
thermodynamic model. 


Martini[23]classified the Stirling analytical model into 
zero-order, first-order, second-order, third-order and fourth- 
order. The zero-order analysis employs empirical relation- 
ships derived from fitting a substantial volume of experimen- 
tal data, facilitating quick and straightforward calculations 
of the Stirling engine’s power output and efficiency. This 
methodology, represented by the Beale number method[24], 
can be beneficial for the engine’s qualitative examination. 
First-order analysis of the Stirling cycle incorporates an ex- 
tensively idealized model, utilizing variables like the temper- 
atures of the cold and hot sources, pressure, rotational veloc- 
ity, piston diameter, and piston stroke. It predicts the out- 
put work and cycle efficiency of the Stirling engine based 
on the conservation of energy between the working spaces. 
Schmidt’s analysis method [24] is a typical example of a first- 
order analysis method, where the overall engine is integrated 
into three spatial segments (dead space, expansion space, and 
compression space). Herein, the prevailing temperatures are 
presupposed to remain constant, often referred to as isother- 
mal analysis. However, the exceedingly ideal nature of this 
analysis results in calculations that deviate significantly from 
real-world Stirling engines. However, the exceedingly ideal 
nature of this analysis results in calculations that deviate sig- 
nificantly from real-world Stirling engines. In contrast, the 
second-order Stirling analysis model elaborates on potential 
power losses and heat dissipations, rendering more precise re- 
sults and gaining widespread usage for Stirling cycle period 
analyses. Nevertheless, there is still a gap between the calcu- 
lation accuracy of traditional second-order analysis methods 
and experimental values, with an absolute error of about 30- 
40% in the estimation of cycle efficiency[24]. To refine pre- 
diction accuracy, numerous studies have aimed at optimizing 
the second-order analysis model. Sayyaadi’s team combined 
finite speed thermodynamics with multiple losses and pro- 
posed Simple II[25] and CAFS[26] models, which reduced 
the error between the simulated and experimental values of 
cycle efficiency to about 6.1%-18.7%. Based on the ideal 
adiabatic model, Ni et al.[27] considered heat losses such as 
incomplete regeneration losses in the regenerator and power 
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Fig. 2. Nuclear fission system and thermal management system of 
Kilopower[34] 


losses such as flow resistance losses and developed the ISAM 
model. The error between the simulated and experimental 
values of cycle efficiency was reduced to about 9.8%-19.9%. 
The third-order Stirling analysis model[28—30] divides Stir- 
ling into multiple nodes in the mainstream direction, solves 
the partial differential equations of mass, momentum, and 
energy conservation at each node to obtain heat and output 
power, and has higher computational accuracy. However, its 
calculation speed is much slower than that of the second-order 
model. The fourth-order analysis method, also known as the 
Computational Fluid Dynamics (CFD) method[31-—33], has 
higher requirements for computing resources and is difficult 
to carry out system analysis. Bearing in mind the need for 
precision and calculation speed to facilitate the subsequent 
development of the thermoelectric conversion module of the 
heat pipe-cooled reactor system analysis code, this paper se- 
lects the second-order analysis methodology for the progres- 
sion of the Stirling analysis model. 

This paper proposes to refine the Simple analysis method 
by incorporating additional loss terms in the loss module, as 
seen in models such as Simple-II, CAFS, and ISAM. These 
additional losses include conductive heat loss, shuttle loss of 
the displacer piston, leakage loss, finite speed loss of the pis- 
ton, and mechanical friction loss. 

It is worth noting that, at present, the free-piston Stirling 
engine is the optimal choice for kilowatt-scale space nuclear 
power systems. In contrast, the application of the GPU-3 in 
this field appears somewhat out of place. However, most 
scholars and research institutions use the GPU-3, a well- 
known beta-type Stirling engine, as a research case for aca- 
demic studies. As early as the last century, NASA conducted 
extensive testing on it, confirming its efficiency and reliabil- 
ity. Notably, Yang et al.[15] proposed a dynamic model of 
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Fig. 3. GPU-3 Stiring engine 


a lunar surface nuclear power system that combines a sec- 
ondary Stirling cycle analysis method, using the GPU-3 Stir- 
ling device. They analyzed the transient response of the re- 
actor under different operating conditions and various space 
environmental temperatures. 


The accuracy of the model in this paper is validated against 
various second-order analysis models and experimental data 
from GPU-3 (Fig.3) and the free-piston Stirling engine RE- 
1000 (Fig.4) [35, 36] . The calculation has proven that the 
model is suitable for simulating both beta-type and free- 
piston Stirling engines. The calculation of the average pres- 
sure incorporates a pressure convergence iterative method, 
enabling rapid convergence of the average pressure. Based 
on the developed model, an appraisal of the operational char- 
acteristics indwelling the Stirling engines is undertaken. This 
paper proceeds to simulate Stirling cycles with usage of hy- 
drogen and helium as working mediums under varying pres- 
sures and rotational velocities to interrogate how their distinct 
thermophysical properties bear upon the engine performance. 
Moreover, this paper probes the effect of regenerator param- 
eters such as porosity, mesh diameter, and the axial length of 
the regenerator on output work and cycle efficiency at differ- 
ent rotational speeds and suggests optimization directions for 
these parameters based on the simulation outcomes. Given 
that the heater and cooler of a Stirling engine are in direct 
contact with the heat source and the cold sink, any changes in 
the temperatures of these sources will directly alter the tem- 
peratures of the heater and cooler, thus affecting the engine’s 
performance. To scrutinize the consequences of temperature 
differentials at the hot and cold ends, Stirling cycles are sim- 
ulated under various temperature ratios. 
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Fig. 4. RE-1000 Stiring engine 


I. STIRLING MODEL 
A. Adiabatic modeling 


Models based on the oscillatory flow characteristics with- 
out pressure gradients are widely used for Stirling cycle char- 
acterization. The Stirling engine is divided into five chambers 
(shown in Fig.5), namely, expansion space, heater, regener- 
ator, cooler, and compression space, and the second-order 
analysis method is based on the ideal cycle, where the mass 
and energy conservation equations and the gas equation of 
state of each chamber are solved to obtain an analytical so- 
lution, which allows further calculation of the various power 
losses and heat losses, and obtains the output power and the 
required heat that are closer to the actual situation. During 
the development of the ideal adiabatic model, the following 
assumptions were made: 

(1) The working fluid is an ideal gas; 

(2) There is no leakage of the working fluid in the cycle, 
and its mass remains constant; 

(3) The temperature of the working gas in the cooler and 
heater is considered to be equal to the wall temperature and is 
held constant; 

(4) The instantaneous pressure throughout the system is 
equal; 

(5) Heat leakage between the compression space and the 
expansion space, as well as heat transfer to the environment, 
is negligible; 

(6) The kinetic and potential energy of the gas flow are 
ignored. 

The adiabatic model differential equation set is shown in 
Table 1, and the subscripts c, k, r, h, and e represent the com- 
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Fig. 5. The diagrammatic representation of the model in the adiabatic analysis methodology 


pression cavity, cooler, return heaters, heaters, and expansion 
cavity, respectively. The subscripts ck, kr, rh, and he repre- 
sent the four interfaces between the five components of the 
compression space, cooler, heat return, heater, and expansion 
space, respectively. 


B. Modification of the Adiabatic Model 


The Simple analysis method only considers the effect of 
non-ideal heat transfer and pressure loss in the heat ex- 
changer. In this paper, based on the Simple analysis method, 
the mass leakage loss term and the distribution piston shuttle 
loss are coupled in the energy conservation equation. 

(1)Shuttle heat loss. Shuttle loss is primarily caused by 
the reciprocating movement of the displacer piston between 
the hot and cold cylinders. There is a significant temperature 
gradient at both ends of the displacer piston, leading to heat 
loss as some of the heat is directly conducted to the cold end 
through the piston. The differential expression for this is: 


2 
dQ shu = IL, (Te a Te) (1) 


where, S is the displacer stroke; k, is the gas thermal conduc- 
tivity; Da is the displacer diameter; J is the gap between the 
displacer and the cylinder. 


I. CALCULATION AND VALIDATION 
A. Calculation process 


(2)Seal leakage loss. In the actual operation of the engine, 
there will be a certain pressure difference between the com- 
pression space and the buffer space, which will lead to a por- 
tion of the gas through the gap between the piston and the 
cylinder wall in the compression space and the buffer space 
back and forth, resulting in leakage losses. The leakage to the 
buffer space was calculated as follows: 
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Table 1. Ordinary differential equations in the Adiabatic analysis 
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Fig. 7. The geometric relationship of the drive equation 


where, P is the pressure in compression space; D is the diam- 
eter of the cylinder; Ppuf fer is the pressure of buffer space; 
J is the gap between piston and cylinder; L, is the length of 
the piston; 


Coupling the shuttle loss and leakage loss with the differ- 
ential equations in the adiabatic analysis method, the pressure 


and pressure differential terms are changed to: 
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B. Non-idea heat transfer 


The heat recovery performance of a Stirling engine is eval- 
uated by the efficiency of the regenerators. The efficiency of 
the regenerator is defined as the ratio of the actual enthalpy 
change to the maximum enthalpy change of the working fluid 
in the regenerators. When the gas flows from the cooler to the 
heater, the temperature at which it exits the return heaters will 
be slightly lower than the temperature of the heater. The gas 
will need to absorb more heat from the heater (external heat 
source), resulting in a less efficient heat cycle. So the actual 
heat absorption and heat release can be written as: 


Qn = Qni + Qrloss = Qni Tr Qri (1 = £) 
Qk am Qki = Qrioss = Qh T Qri (1 = E) 


where, Qn is heat addition to the working fluid in the heater; 
Qk is heat rejection to the cooler; Qross is heat loss in the 
imperfect regenerator; Qr; is the amount of regenerator heat 
transferred in ideal process and e is the effectiveness of the 
regenerator calculated as: 
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where, NTU is the number of heat transfer units; Awg inter- 
nal wetted area of the regenerator; A is cross-sectional area 
of the regenerator and St is the Stanton number calculated as: 


St = 0.46Re~°* Pr ®t (9) 


where, Re is the Reynolds number; Pr is the Prandtl num- 
ber.In this paper; Pr is taken as a constant of 0.71. The 
Reynolds number in the regenerator is calculated as: 
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l-9 
where, dm is the diameter of the regenerator wire; y is the 
porosity coefficient. The heater and cooler are also non-ideal, 
and the actual temperature of the gases in the two heat ex- 
changers depends on the heat losses in the return heaters and 
the convective heat transfer between the gases and the walls 


in the heater and cooler. The actual gas temperature in the 
heat exchanger is given as follows: 
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where Qacn is the actual heat absorbed in the heater; Qack 
is the actual heat rejection of the cooler; Wac is the actual 
indicated power output; h is the heat transfer coefficient in 
the heater and cooler. The heat transfer coefficient is given as 
follows: 

fruc 

h= P 13 

2dPr (13) 
In the analysis of the Stirling engine cycle, Twn and Tug de- 
note the wall temperatures of the heater and the cooler, re- 
spectively, and are generally set as constant values. 


C. Other losses effects 


(1)Thermal conductivity loss. The regenerator is connected 
to the cooler and the heater, and the heat from the hot end 
flows through the cylinder wall of the return heaters to the 
cold end, which is the main source of heat conduction loss. It 
can be evaluated according to the Fourier law[37]: 
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(14) 
where, km is the heat conduction coefficient of the material; 
Aw is the Thermal conductivity cross section; / is the conduc- 
tivity length. 

(2)Flow resistance loss. The pressure loss is caused by fric- 
tion due to the viscosity of the fluid, and the loss caused by 
pressure drop is known as flow resistance loss. The multi- 
layered annular mesh structure within the regenerator is the 
primary cause of flow resistance loss. Pressure drop and flow 
resistance loss can be calculated using semi-empirical formu- 
las: 
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where, Cref is the Reynolds friction factor; ju is the dy- 
namic viscosity; V is the volume of the working area; d is 
the hydraulic diameter; A is the flow cross-section area.The 
Reynolds friction factor C’-¢¢ is caculated as: 
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(3)Finite speed of piston and mechanical friction loss. Rel- 
evant studies have indicated that the instantaneous pressure 
on the piston surface differs from the instantaneous average 
effective pressure within the engine cylinder[38]. The actual 
work of compression and expansion is different from the the- 
oretical calculations based on classical thermodynamics. This 
pressure loss is caused by the finite speed motion of pressure 
waves generated by the piston within the working space. The 
formula for calculating the work loss due to finite speed of the 
piston and mechanical friction can be calculated as follows: 
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where, Pm is the mean effective working pressure; up is the 
linear speed of the piston; f is the frequency; AP, is the 
pressure loss caused by mechanical friction of the compo- 
nents. The + represents the compression process (+) and the 
expansion process (-), respectively. Pm can be calculated as : 
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where, rẹ, is the compression ratio, ry, = Pae, 


The algorithm for the model of this paper is illustrated in 
Fig.6. In the first step, initialize Stirling engine geometry, 
hydrodynamic and thermodynamic parameters. Geometric 
equations for V., Ve, dV. and dV. as a function of the an- 
gle of rotation are obtained from the configuration of the en- 
gine. Calculate and initialize hydrodynamic parameters such 
as cross-section area, wet area, hydraulic diameter, etc. based 
on geometric parameters. Initialize thermodynamic parame- 
ters such as specific heat capacity, gas constant and thermal 
conductivity; 

Secondly, initialize the computational parameters. The 
temperatures of the expansion and compression space are ini- 
tialized to the wall temperature Twp and T;,,,. The initial mass 
is computed by means of a first-order model; 

Thirdly, initialize angle @ as 0° and set the calculation step 
to 1°. The angle is initialized at the end of each cycle, and 
for the initial value problem, the system of differential equa- 
tions is solved using the fourth-order Lunge-Kutta method 
with 9=360° to complete the calculation of one cycle, and 
it is judged whether the cavity temperature at the end of the 
cycle is equal to the temperature at the beginning of the cycle 
until the gas temperature reaches the limit of convergence to 
turn on the next step of the calculation; 

Then, calculate the regenerator’s regeneration efficiency 
and the convection heat transfer coefficients of the heater and 
cooler, for details refer to literature[24]. Based on the ideal 
heat exchange amount obtained from the previous step, derive 
the actual heat exchange quantity between the heat exchanger 
wall and the working fluid. Compute the new temperature of 
the heat exchanger and use it as the initial value to return to 


the third step of the ideal adiabatic model for calculation. Re- 
peat this process until the convergence criteria are met, and 
then proceed to the fifth step of the calculation. 

Finally, calculate the other losses to obtain the actual heat 
and output work. The method for solving the classical fourth- 
order Runge-Kutta method is shown as follows. 
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D. Average pressure convergence 


Here, an iterative method for pressure convergence is ap- 
plied. During the process of solving the Stirling cycle, it is 
necessary to calculate the average pressure to meet the condi- 
tion of average pressure convergence. The average pressure 
in the compression space over a cycle should be equal to the 
set average pressure, that is: 
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perepresents the residual value for determining pressure con- 
vergence. After the first cycle of the solution is completed, 
the initial pressure of the model is adjusted until the calcu- 
lated average pressure meets the periodic stable solution as 
described by the above formula, that is: 


Pstart_new = Pstart_last — PePstart_last (22) 


Employing this iterative method allows for the rapid iden- 
tification of a solution that meets the target average pressure 
across the entire system. 


E. Verification and validation 
1. Validation of the GPU-3 Engine 


The GPU-3[23, 39] Stirling engine, a portable generator 
set developed by General Motors for the Army, has well- 
established equations of motion for the rhombic drive mech- 
anism to calculate the expansion and compression space vol- 
ume changes, as shown in Fig.7. The detailed design dimen- 
sions and parameters for the GPU-3 are presented in Table 
2: 

To validate the model, simulations (Helium, Twh = 922 
K, Twk = 286 K) were carried out against the GPU-3 Stir- 
ling engine experiment with two different pressure conditions 
(Pmean = 2.76 MPa and Pmean = 4.14 MPa). Meanwhile, the 
computational results are compared with the simulation data 
of available second-order analysis models (Simple-II, CAFS, 
ISAM, Simple), the third-order analysis software Sage. The 
experimental data were obtained from NASA experiments on 
the GPU-3 Stirling engine. 


Table 2. GPU-3 engine dimensions and parameters. 


Parameters Values Parameters Values 

Clearance volumes Swept volumes 
Compression space 28.68 cm? Compression space 113.14cm? 
Expansion space 30.52 cm? Expansion space 120.82 cm? 

Heater Cooler 
Tube number 40 Tube number 312 
Tube inside diameter 3.02 mm Tube inside diameter 1.08 mm 
Tube length 245.3 mm Tube length 46.1 mm 
Void volume 70.38 cm? Void volume 13.8 cm? 

Regenerator Drive 
Void volume 50.55 cm? Connecting rod length 46 mm 
Length 22.6 mm Crank radius 13.8 mm 
Internal diameter 22.6 mm Eccentricity 20.8 mm 
No. per cylinder 8 Displacer stroke 31.2 mm 
Diameter of wire 0.04 mm Internal diameter of cylinder 69.9 mm 
Porosity 0.697 Working fluid Helium 
Material Stainless steel Frequency 41.72 Hz 

Table 3. RE-1000 engine dimensions and parameters. 

Parameters Values Parameters Values 

Heater Cooler 
Tube number 34 Fin number 115 
Void volume 27.33 cm? Void volume 20.43 cm? 

Regenerator Volume 
Porosity 0.759 Expansion space 27.74cm? 
Diameter of wire 0.0889 mm Compression space 54.80 cm? 
Void volume 56.37 cm? 

Displacer Phase angle S15” 
Stroke 24.5 mm Frequency 30 Hz 
Diameter 56.4 mm Working fluid Helium 

Power piston Heater temperature 873 K 
Stroke 28 mm Cooler temperature 297 K 
Diameter 57.2 mm Mean pressure 70.6 MPa 
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Fig.8 and Fig.9 illustrates the relationship between the 


working fluid temperature variation with angle in the expan- 
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Fig. 9. Variation of the working fluid temperature 


sion and compression cylinders, as well as the P-V (pressure- 
volume) diagram of the cycle. The temperature of the cavity 
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Fig. 10. Comparison of output power between different models and experimental values (Helium) 
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Fig. 11. Comparison of thermal efficiency between different models and experimental values (Helium) 


returns to the initial point after one complete rotation of the 
Stirling rhombic drive mechanism, and the cycle repeats peri- 
odically, which is in accordance with the basic characteristics 
of the Stirling cycle. 

Fig. 10 presents a comparison of the output work results be- 
tween the model in this paper and the existing second-order 
model at various frequencies. Fig.11 displays a comparison 
of the cycle efficiency results between the model in this paper 
and the existing second-order model across different frequen- 
cies. Compared to the Simple model (Pmean=4.14 MPa), the 
maximum power error of this paper’s model has been reduced 
by 152.5%, and the maximum efficiency error (as difference) 
has been decreased by 10.67%. The average reduction in 
power error across all frequency conditions is 110.0%, and 
the average reduction in efficiency error is 8.43%. Fig.12 
illustrates the comparison of simulation values between the 
model in this paper and the existing second-order model un- 
der the standard operating conditions of 4.14 MPa pressure 


and 41.67 Hz frequency. In comparison to the Simple model, 
the simulation error for output work has been lowered by 
46.2%, and the simulation error (as difference) for cycle ef- 
ficiency has been reduced by 8.7%. In comparison with the 
third-order analytical model Sage, the output work exhibited 
maximum and minimum errors(as difference) of 28.3% and 
4.5%, respectively, with an average error of 11.5%. The cy- 
cle efficiency demonstrated maximum and minimum errors 
(as difference) of 5.43% and 0.30%, respectively, averaging 
an error of 1.55%. 

Comparing the simulated and experimental values allows 
for the analysis of some working characteristics of the Stir- 
ling engine. When helium is used as the working fluid, the 
output work of the Stirling engine decreases as the frequency 
increases under high-frequency conditions (with a working 
frequency greater than 41.72 Hz), and the cycle efficiency 
shows an overall downward trend within the frequency anal- 
ysis range we have defined, with this downward trend being 
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Fig. 12. Bar Chart of Simulation Results Comparison. 


Table 4. Comparison of Experimental and Simulated Values of the RE-1000 engine. 


Parameters Present study Experiment[40] 
Heat in(W) 3940 (-2.4%) 4038 
Indicated power(W) 1138 (10.4%) 1030 
Indicated efficiency(%) 28.9 (4.4%) 25.5 


particularly pronounced in the high-frequency region. On the 
contrary, in the subsequent parameter sensitivity analysis, it 
was found that when hydrogen is used as the working fluid, 
both the output work and cycle efficiency of the Stirling en- 
gine exhibit a significant upward trend. This difference is 
caused by the different material properties of the two sub- 
stances. Additionally, as the working frequency increases, 
the error between the experimental and simulated values in 
the engine also gradually increases, which is due to other un- 
determined losses caused by high-frequency oscillating flows. 


2. Validation of the RE-1000 Engine 


In order to further validate the correctness of the model, 
experimental data from the RE-1000 free-piston Stirling en- 
gine, developed by NASA, was utilized for the validation pro- 
cess. The model was validated against the operating condi- 
tions and experimental outcomes specified in reference [40] 
for test case #1011. Several studies have referenced this data 
for the validation of the RE-1000 engine [35, 36]. The de- 
tailed design dimentions and parameters of the RE-1000 are 
presented in Table 3. Due to the periodic motion of the dis- 
placer and the power piston, the volume changes of the com- 
pression space and the expansion space are described as fol- 
lows: 


Ve = Veo = AgXa sin (wt) 
Ve = Veo + (Aa — Aroa) Xa sin (wt) (23) 
— A,X> sin (wt — a) 


where, the terms Ve and V, denote the initial volumes of 
the expansion and compression spaces, respectively, at the 
equilibrium positions of the displacer and the power piston; 
Aq, Ap and Aroq represent the cross-sectional areas of the 
displacer, power piston, and connecting rod, respectively; Xq 
and X, are the displacer amplitude and power piston ampli- 
tude; a is the angular difference between the motions of the 
displacer piston and the power piston. It is used to describe 
the synchronization between the pistons’ movements and can 
impact the efficiency and power output of the engine. 


Table 4 shows the results of the comparison between exper- 
imental and simulated values. The relative errors for the input 
heat, output work, and cycle efficiency are correspondingly - 
2.4%, 10.4%, and 4.4%. The validation against experimental 
and simulated data from the GPU-3 and RE-1000 Stirling en- 
gines has demonstrated the model’s high level of precision. 
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Fig. 13. The variation of output power with frequency (Hydrogen) 
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Fig. 14. The variation of efficiency with frequency (Hydrogen) 


IV. PARAMETER SENSITIVITY ANALYSIS 
A. Effect of frequency and different working fluids 


The selection of working fluids for Stirling engines is pre- 
dominantly focused on gases, with commonly used working 
fluids being Nə, Hə, CO2, He, and air. Due to the supe- 
rior heat exchange characteristics of small molecular weight 
gases, most technologically advanced high-power Stirling en- 
gines employ He or H3 as their working fluids. 


In order to investigate the effects of frequency and dif- 
ferent work materials on the output work and cycle effi- 
ciency, hydrogen is selected as the work material (Twh=977 
K, Twk=288 K), with two different pressure variables set 
at 1.38 MPa and 2.76 MPa. Fig.l3a and Fig.13b give the 
changes of the output work with the frequency. Fig.14a and 
Fig.14b give the changes of the cycle efficiency with the fre- 
quency. 
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Fig. 15. Effect of rotation speed on Reynolds number of regenerator 
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Fig. 16. Effect of rotation speed on pressure drop of regenerator 
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Fig. 17. The maximum pressure drop of the heat exchanger varies 
with frequency 


In scenarios of slower rotation, the efficiency of the sys- 
tem increases with increased rotation speed. However, as the 
speed of rotation continues to rise, there is a decrease in sys- 
tem efficiency, resulting in a general trend of initial increase 
followed by a decrease. When the work material is H2, there 
is a clear rise in the power output of the Stirling engine with 
the increase in rotational speed, and the upward trend is sharp. 
On the other hand, when the work material is He, the faster 
rotation speed clearly has a detrimental impact on the Stir- 
ling engine’s output, with a less steep upward trend. Given 
that the regenerator embodies a screen-like structure, a sub- 
stantial pressure drop transpires as the working mass flows 
through it, consequently yielding a loss attributed to flow re- 
sistance. Fig.15 and Fig.15 demonstrates the effect of differ- 
ent rotation speeds on the pressure loss and Reynolds num- 
ber of regenerator. Fig.17 demonstrates the variation of the 
maximum pressure drop within the heat exchanger with fre- 
quency. With the increase in frequency, the Reynolds num- 
ber and pressure loss inside the heat exchanger both increase. 
Compared to the pressure loss in the regenerator, the pres- 
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sure losses in the cooler and heater are negligible, confirming 
that the wire mesh structure of the regenerator is the primary 
source of pressure loss. The high rotational speed of the Stir- 
ling engine makes the oscillating flow and heat exchange of 
the work material inside the body tend to be complicated, re- 
sulting in a significant increase in the flow resistance loss or 
local loss, which in turn affects the cycle efficiency. Hə pos- 
sessing a lower dynamic viscosity coefficient, produces rel- 
atively smaller flow resistance loss when operating at high 
rotational speeds than He. However, Ho exhibits poor seal- 
ing properties, making it prone to leakage and explosion, and 
it can cause issues such as hydrogen embrittlement in cer- 
tain materials. Therefore, helium is commonly chosen as the 
working fluid. 


B. Effect of regenerator parameters 


The regenerator features an internal structure filled with a 
metal wire mesh. Porosity is represented by the ratio of the 
regenerator’s porous volume to its overall volume. Fig.18 de- 
picts the variation of output work and cyclic efficiency with 
changes in porosity and wire mesh diameter. For metal wire 
meshes of different diameters, the cyclic efficiency tends to 
increase initially and then decrease. This trend is attributed to 
the fact that an appropriate increase in porosity enlarges the 
hydraulic diameter and reduces the wetted area of the mesh 
in contact with the working fluid. This reduction in pressure 
drop across the regenerator leads to a decrease in flow resis- 
tance losses, thereby enhancing the output work and cyclic 
efficiency of the engine. At lower porosities, a denser mesh 
configuration results in a smaller hydraulic diameter and an 
increased wetted area for the regenerator, which can signif- 
icantly impact the flow resistance and adversely affect the 
cyclic efficiency and output work. 

When the porosity is below 0.797, Stirling engines 
equipped with larger diameter wire meshes have higher cyclic 
efficiencies compared to those with smaller diameter meshes. 
This is due to the larger convective heat transfer area provided 
by the larger diameter meshes under conditions of low poros- 
ity. On the contrary, when the porosity surpasses 0.797, wire 
meshes of a smaller diameter demonstrate superior cyclic effi- 
ciencies, attributable to their comparatively lesser flow resis- 
tance losses. At a wire mesh diameter of 0.04 mm, a turning 
point in cyclic efficiency is observed at a porosity of 0.797. 
With the further escalation of porosity to 0.847, there is a 
slight rise of 0.1% in output power, while the cyclic efficiency 
drops by 3%. In conclusion, the engine’s output work and 
cyclic efficiency are most optimal at a porosity of 0.797. 

Fig.19 illustrates the variation of output work and cycle ef- 
ficiency with the length of the regenerator. It can be observed 
from the figure that at a mean pressure of 4.14 MPa, an in- 
flection point appears in the cycle efficiency curve when the 
regenerator length is around 0.025 meters. When the regen- 
erator length is less than 0.025 meters, an increase in length 
leads to a larger heat exchange area, allowing for more suffi- 
cient heat transfer within the regenerator’s metal wire mesh, 
which in turn increases cycle efficiency with the increase in 
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Fig. 18. The influence of the porosity and wire diameter 
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Fig. 19. The influence of the regenerator length 


length. Moreover, since the pressure losses due to flow re- 
sistance will inevitably increase with length, the output work 
will necessarily decrease as the regenerator lengthens. This 
creates a situation where, when the regenerator length ex- 
ceeds 0.025 meters, the positive effect of increased heat ex- 
change area on cycle efficiency is not sufficient to offset the 
negative effect caused by increased losses, resulting in a de- 
cline in cycle efficiency. It can be observed that under the con- 
dition of Pmean=4.14 MPa, the optimal regenerator length is 
around 0.025 m. 


C. Effect of temperature of heat exchangers 


In the actual operation of a space nuclear reactor, the hot- 
end temperature of a Stirling system fluctuate with alterations 
in the core temperature of the nuclear reactor, and the per- 
formance of the radiative heat rejection system can also af- 
fect the cold-end temperature. Shifts in both the hot-end and 
cold-end temperatures directly impact the quantity of heat 
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Fig. 20. The influence of the temperature ratio 


absorbed and released by the working fluid gas, thereby in- 
fluencing the output work and cyclic efficiency of the Stir- 
ling engine. Consequently, with respect to the space nuclear 
power reactor, it is necessary to assess the impact of cold-end 
and hot-end temperatures on the performance of the Stirling 
system. 


Fig.20 demonstrates the variation in output work and cyclic 
efficiency of the Stirling engine at different temperature ratio 
conditions as the cooler temperature increases from 280 K 
by increments of 20 K up to 440 K. At a temperature ratio 
of 0.25, the output work increases from 4874 W to 5301 W, 
marking an increment of 2.3%. At a temperature ratio of 0.35, 
the output work rises from 3146 W to 3480 W, an increase 
of 5.4%. When the temperature ratio is 0.45, the output work 
surges from 1773 W to 2037 W, denoting an increase of 9.8%. 
When the temperature ratio is fixed, alterations in the cooler 
temperature negligibly influence the output work and cyclic 
efficiency. In actual scenarios, due to material constraints, 
the heater of a Stirling engine cannot reach excessively high 
temperatures, so the ideal case is considered here. With a 


constant cooler temperature, when the temperature ratio in- 
creases from 0.25 to 0.35, the output work decreases on aver- 
age by 34.8%, and the cyclic efficiency decreases by 8.61% 
on average. When the temperature ratio increases from 0.35 
to 0.45, the output work decreases on average by 42.3%, and 
the cyclic efficiency decreases by an average of 10.0%. It can 
be observed that the lower the temperature ratio, the higher 
the system’s output power and cyclic efficiency. Conversely, 
the higher the temperature ratio, the lower the system’s output 
power and cyclic efficiency. With a constant cooler tempera- 
ture, a higher heater temperature is more conducive to achiev- 
ing satisfactory output work and cyclic efficiency. Similarly, 
under the condition of constant heater temperature, the lower 
the cooler temperature, the better the engine performance. 
However, in practice, due to the impact of environmental fac- 
tors, the cooler cannot achieve excessively low temperatures, 
and this should be taken into account based on actual condi- 
tions. The operating temperature of the working fluid in high- 
temperature heat pipes exceeds 730K. Among them, the most 
frequently used sodium heat pipes have a working tempera- 
ture range of 700-1100K. The hot end of the Stirling engine 
transfers heat to the high-temperature heat pipe through solid- 
to-solid contact, while the cold end can only dissipate heat 
into the space environment through radiative cooling. There- 
fore, for space power devices, thermal efficiency is particu- 
larly important; a higher thermal efficiency means that under 
the same weight and cost conditions, it allows for a higher 
power output or a longer operational time. 

In the parameter sensitivity analysis of the cold and hot 
end temperatures, we have fixed the temperature boundaries 
to simulate the operating conditions under steady-state. How- 
ever, in practical engineering applications, temperature con- 
trol and thermal management are typically complex and crit- 
ical issues. For example, during the reactor start-up process, 
there is a significant concern about controlling the thermal 
balance throughout the entire process, from the heat gener- 
ated in the core to its transfer via high-temperature heat pipes 
to the Stirling engine and finally to the heat rejection system. 
Heat pipes are often quite fragile, and initiating a cold start 
of the Stirling engine could lead to thermal imbalance, caus- 
ing the heat pipes to fail and posing a risk to the safety of 
the reactor. Therefore, for the start-up issue, it is necessary 
to devise a reasonable temperature control plan to balance the 
heat transfer between various components. We will explore 
this issue in subsequent work. 


V. CONCLUSION 


Stirling engines, with their high efficiency, strong reliabil- 
ity, and compact structure, have attracted significant atten- 
tion from researchers, especially in recent decades, leading 
to substantial development and an irreplaceable role in var- 
ious fields. To investigate the operational characteristics of 
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Stirling engines, this paper establishes a second-order anal- 
ysis model, thoroughly validates it, and conducts extensive 
parameter analysis based on the validated model. The spe- 
cific work includes: 

(1) The model developed in this paper is fully validated 
against the experimental data of the GPU-3 Stirling engine 
from NASA. Under the condition of 4.14 MPa pressure, com- 
pared to the Simple model, the maximum power error has 
been reduced by 152.5%, and the maximum efficiency error 
has been lowered by 10.67%. The average reduction in power 
error across all frequency conditions is 110.0%, and the av- 
erage reduction in efficiency error is 8.43%. Furthermore, 
comparisons were also conducted with the experimental data 
of the RE-1000 free-piston Stirling engine. The results show 
that the model developed in this paper demonstrates equiv- 
alent accuracy compared with existing second-order models, 
thereby confirming the correctness of the model. 

(2) This study investigates the impact of working fluids on 
the operational characteristics of Stirling engines. The results 
indicate that hydrogen, with its low dynamic viscosity coeffi- 
cient and reduced flow resistance losses, can generate higher 
output power when used as the working fluid in Stirling en- 
gines. Additionally, the research reveals that high rotational 
speeds lead to complex oscillatory flow phenomena, resulting 
in increased pressure drops. Furthermore, as the Reynolds 
number increases, flow resistance losses rise, leading to a re- 
duction in the engine’s output power. 

(3) The output power and cycle efficiency of the Stirling cy- 
cle were simulated under different wire mesh porosities and 
diameters. The results showed porosity, wire diameter, and 
regenerator length all have a significant impact on the perfor- 
mance of the Stirling engine. When the wire mesh diameter 
is 0.04 mm, the optimal porosity is 0.797, and the optimal re- 
generator length is around 0.025 m. 

(4) The performance of a Stirling engine in a space nuclear 
reactor is significantly influenced by the hot-end and cold- 
end temperatures. Variations in these temperatures directly 
affect the heat exchange and efficiency of the engine. Opti- 
mal temperature ratios enhance output work and efficiency, 
while extremes can lead to significant reductions. Thermal 
management, particularly during start-up, is critical to pre- 
vent thermal imbalance and ensure reactor safety. . 
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